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Abstract 

We investigate the marginally stable modes of the scalar (vector) pertur- 
bations in the AdS soliton background coupled to electric field. In the probe 
limit, we find that the marginally stable modes can reveal the onset of the 
phase transitions of this model. The critical chemical potentials obtained from 
this approach are in good agreement with the previous numerical or analytical 
results. 
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1 Introduction 



The AdS/CFT correspondence [1] provides a powerful theoretical method to understand 
the strongly coupled field theories. Recently it has been applied to study the holographic 
models of superconductors (superfluids) phase transitions since the work [2, 3]. For reviews 
see [4] . Since the condensed matter physics deals with the systems at finite charge and finite 
temperature, from the AdS/CFT correspondence the dual gravity should be described by 
a charged black hole. 

The phase transition studied in [2, 3] is actually a holographic superconductor/metal 
phase transition. The model for holographic superconductors can be simply constructed 
by an Einstein-Maxwell theory coupled to a complex scalar field. In particular, when 
the temperature of the black hole is below a critical temperature, the black hole solution 
becomes unstable to develop a scalar hair near the horizon. And the condensation of the 
scalar hair breaks the U(l) symmetry of the system. From the AdS/CFT correspondence, 
the complex scalar field is dual to a charged operator in the boundary field theory. And 
the breaking of the U(l) symmetry in gravity causes a global U(l) symmetry breaking in 
the dual boundary theory. This induces a superconductor (superfluid) phase transition [5]. 

The holographic insulator /superconductor phase transition was first studied in [6]. In 
particular, they used a five-dimensional AdS soliton background [7] coupled to a Maxwell 
and scalar field to model the holographic insulator/superconductor phase transition at 
zero temperature. The normal phase in the AdS soliton is dual to a confined gauge theory 
with a mass gap which resembles an insulator phase [8]. When the chemical potential is 
sufficiently large beyond a critical value, the AdS soliton becomes unstable to form scalar 
hair which is dual to a superconducting phase in the boundary field theory. The holographic 
insulator /superconductor phase transition was also studied in [9, 10, 11, 12, 13]. 

To reveal the stability of a spacetime background, a powerful method is to study the 
quasinormal modes (QNMs) of the perturbations in this background, for reviews see [14, 
15, 16]. If the imaginary part of the QNMs is negative, the mode will decrease in time 
and the perturbation will finally disappear which indicates that the background is stable 
against this perturbation. On the contrary, if the imaginary part of the QNMs is positive, 
this implies that the background is unstable against this perturbation. The interesting 
thing is that if the perturbation has a marginally stable mode, i.e. u — 0, one always 
expects that this is a signal of instability, or rather a phase transition may occur. For the 
detailed discussions of marginally stable modes, one can refer to [2]. 

In this paper, we employ the idea of marginally stable modes to study the s-wave 
and p-wave holographic insulator/superconductor phase transitions in AdS soliton back- 
ground at zero temperature. The construction of the system is like the one in [6]. We 
take advantage of Horowitz and Hubeny's method [17] to study the QNMs of the scalar 
(vector) perturbations of this system. By increasing the chemical potential from zero to 
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some critical value, the marginally stable modes will turn out. This means at the crit- 
ical chemical potential the AdS soliton background becomes unstable and will prefer to 
be an AdS soliton background coupled with nonzero charged scalar (vector) fields. This 
argument is consistent with the previous studies on holographic insulator / superconductor 
phase transitions [6, 9, 10]. In particular, the critical chemical potentials we get from the 
marginally stable modes are in good agreement with the ones in [6, 9, 10]. Actually, there 
are multiple marginally stable modes corresponding to various critical chemical potentials. 
These modes are related to the modes of node n = 1, 2, 3 • • ■ . However, they are unstable 
due to the oscillations of scalar (vector) field in the radial direction [2, 18]. By making use 
of the alternative approach, viz. "shooting" method, we plot the behavior of the scalar 
(vector) fields depending on the radial direction. From these diagrams, one can intuitively 
see the "nodes" of these fields. The studies of QNMs in AdS soliton background are also 
investigated in [19, 20]. 

The paper is organized as follows. In Sect. (2), we study the marginally stable modes 
of s-wave field in the probe limit. We find that the critical chemical potentials we derived 
are consistent with the results obtained by previous works. By making use of the same 
procedure, we explore the QNMs of p-wave field in the AdS soliton background in Sect. (3). 
We draw our conclusions and make some discussions in Sect. (4). 



2 S-wave perturbations 

Following Ref.[6], we set up this model in the AdS soliton background [7]: 

dv^ 

ds 2 = L 2 -—- + r 2 (-dt 2 + dz 2 + dy 2 ) + f(r)d X 2 . (1) 
f( r ) 

where, f(r) = r 2 — r^fr 2 and L is the radius of AdS spacetime. In fact, this soliton 
solution can be obtained from a five- dimensional Schwarzschild-AdS black hole by making 
use of two Wick rotations. The asymptotical AdS space-time approaches to a topology of 
R 1 ' 2 x S l near the boundary. And the Scherk-Schwarz circle x ~ X + 7r -^/ r o is required in 
order to have a smooth geometry. The geometry looks like a cigar whose tip is at r = r . 
Because of the compactified direction x, this background provides a gravity description 
of a three-dimensional field theory with a mass gap, which resembles an insulator in the 
condensed matter physics. The temperature in this background is zero. 

It is well known that in this AdS soliton background a simple solution for Maxwell 
gauge field is A t = Const. = fi. Instead of A t = at the horizon required by the AdS black 
holes, A t can be any non-singular value at the tip of the AdS soliton. 

In the probe limit, we introduce a charged scalar field ip as a probe into this background 
which is a neutral AdS soliton with a constant electric potential. The Lagrangian for the 
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charged scalar field is 

^matter = ~| ~ ~ m 2 \lp\ 2 . (2) 

The Euler-Lagrange equations of motion (EoMs) for ip is 

( V M - iqAp) ( V M - iqA^ip - m 2 ip = 0. (3) 

In the following, we assume that ip is real and ijj = F{t,r)H(x)Y(z,y). Substituting tf> 
into Eq.(3) and making the separation of the variables we can reach 

d 2 F(t,r) t ,3 d r f dF(t,r) _ L 2 d 2 F(t,r) 2iq^iL 2 dF(t,r) 



dr 2 r f dr fr 2 dt 2 fr 2 dt 

+^(?V - m 2 r 2 - - e)F(t, r) = 0. (4) 

where A and £ are the eigenvalues of the following equations respectively 

d 2 H( X ) 



dx 

d 2 Y(z,y) , d 2 Y( Zl y) 



. yH( X ) = 0, (5) 
+ eY(z,y) = 0. (6) 



dz 2 dy 2 

where A = 2r n/L, n G Z due to the periodicity of H(x) = H(x + n L/r ) and (GZ. For 
simplicity we will set A = £ = which means that there are no momenta in the (z, y, in- 
directions. 



2.1 Critical behavior from marginally stable modes 



Further, we define F(t,r) = e %u>t R{r) and Eq.(4) becomes (we have set L = 1) 



R"(r) + (^ + -)R'(r) 
J r 



+ + <?/-0 - m 2 r 2 ]R(r) = 0. 



(7) 



where a prime denotes the derivative with respect to r. 

In order to investigate the phase transitions of this model, we recall that the marginally 
stable modes can to some extent reveal this critical behavior [2]. In the studies of QNMs, 
marginally stable modes correspond to oj = which indicates that the phase transition or 
the critical phenomena may occur. We will take advantage of the Horowitz and Hubeny's 
method [17] to study these QNMs. 

It is convenient to convert the r-coordinate to x-coordinate, where x = 1/r. Therefore, 
the infinite boundary is now at x = while the tip is at X — Xq — l/r . In terms of this 
new coordinate x, Eq.(7) becomes 

x 4 d x f(x) 



x d xx R(x) 



— x 



/(*) 



d x R(x) 



/(*) 



x 2 (u 



q/j,) — m 



R(x) = 0. 
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Following the steps of Horowitz and Hubeny [17], we can multiply f(x)/(x — xq) to both 
sides of the above equation, and we reach 



S{x)d 2 x R{x) 

X — Xq 

where the coefficient functions are given by 

f(x)x 4 



' ] -d x R(x)+ . U{X \ 2 R(X) = 1 



(x - x y 



S(x) 
T(x) 
U(x) 



X — Xq 

-fx 3 + x A d x f(x), 
x {lo + q/jL) — m 



(X - Xq). 



(9) 

(10) 

(11) 
(12) 



Note that x = xq is a regular singular point of S(x),T(x) and U(x), and we can polyno- 
mially expand them to a finite order like 



M 



x - x j 



(13) 



n=0 



where M is a finite integer. The series expansion of T(x) and U(x) can be similarly reduced. 

Unlike the ingoing boundary conditions of scalar field near a black hole horizon, the 
boundary conditions here can be a finite quantity at the tip of the AdS soliton. We can 
expand R(x) = (x — Xo) a and substitute it into Eq.(9). Then to the leading order we get 



a(a — l)s + at + u = —Ax^a 2 = 
This corresponds to looking for a solution of the form 



a = 0. 



A' 



R(x) = lim > a 



[X - X ) 



(14) 



(15) 



n=0 



Substituting (15) and (13) into (9) and comparing the coefficients of (x — xo) n for the same 
n we find that 



— 22[k(k - l)s n _ fc + kt n _ k + u n _ k ]a k , 



k=0 



n{n — 1)sq + nto + uq = —AxqU 2 



(16) 
(17) 



We set «o = 1 due to the linearity of Eq.(9). The boundary conditions for the scalar field 
at x = is 



/v 



-R(O) = ^lim a n 



-Xq) 



0. 



n=0 



And the algebraic equation (18) can solve the modes oj. 

In the following numerical calculations, we restrict x$ = 1 and q = 1 just like the ones 
in [6]. In practice, we will expand R(x) to a large order which is N = 300. In order to 
find the marginally stable modes u = of the system, we should restrict our attention to 
the potential U(x) in (12). One finds that u and fi are symmetric (q = 1). This means 
that when fi = the lowest-lying modes of u will exactly be identical to the lowest-lying 
critical chemical potentials \i c when uj = 0. 1 Using this trick, we can easily find the critical 
chemical potentials where the marginally stable modes arise. 




l.o h , , , , i , , , , i , , , , i , , , , d 

-4 -3 -2 -1 

m 2 

Figure 1: The critical chemical potentials for the marginally stable modes versus m 2 of the 
scalar field. The blue curve is taken from the analytical results in [13] while the red points 
are numerically obtained from the QNMs. They are perfectly matched. 

Fig.(l) shows the critical chemical potentials /i c 's for various squared mass of the scalar 
field. The red points are obtained from the numerical calculations of the marginally stable 
modes while the blue curve is taken from the analytical results we have previously derived 
in [13]. They agree with each other perfectly. In particular, when m 2 = —15/4, \x c rs 1.885 
obtained from the marginally stable modes is in good agreement with the result \xi = 1.88 
which is the critical value for the onset of the holographic insulator/superconductor phase 
transitions obtained in [6]. This implies that the marginally stable modes can indeed reveal 
the phase transitions of this model. 

Table (1) shows the first three lowest-lying critical chemical potentials for various m 2, s. 
The /i c 's of the overtone number n = are the critical chemical potentials shown in Fig.(l). 
Other critical chemical potentials of overtone numbers n = 1, 2 can also make the QNMs to 

1 Here, we only take care of the QNMs with positive real part. Therefore, the lowest- lying modes of uj 
represent the modes which have the minimal or less minimal positive real parts. 
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Table 1: The first three lowest-lying critical chemical potentials /i c for various mass squares 
obtained from the calculation of the marginally stable modes. The overtone numbers of 
them are n = 0, 1, 2 from the top to the bottom. 



m 2 = —4 


m 2 = -15/4 


m 2 = —3 


m 2 = -2 


1.373 


1.885 


2.396 


2.815 


3.658 


4.226 


4.792 


5.246 


6.029 


6.603 


7.189 


7.655 



be marginally stable. However, they are expected to be unstable which can be understood 
after the next subsection. The nodes n = 0, 1,2 can also be intuitively seen in the next 
subsection. 

2.2 Critical behavior from the "shooting" method 

Actually there is another way to study the critical behavior of this phase transition. It 
was called the "shooting" method which was commonly used in the previous studies on 
holographic superconductors [4]. Here, we will concisely describe how to make use of this 
"shooting" method to study the critical behavior in AdS soliton background and compare 
it with the results of "marginally stable modes" method. 

We will work in the approximation that A t = Const. = \i and the scalar field only 
depends on r-direction as well as that it is too small to back-react the background. The 
EoMs of i/j(r) is 

f 3 a 2 m 2 

+ (j + p^'( r ) + (£j ~ yMO = 0. (19) 

The boundary condition of ip(r) at the tip is 

rj} = a + b(r - r ) H , (20) 

and near the infinite boundary ip(r) behaves as 

^ r 2-V4+m 2 r 2+V4+TO 2 ^ ' 

In the following calculations, we will set ip^ = in order to turn off the effect of the source 
on the boundary field theory. 2 



2 It is well known that when < a/4 + m? < 1 the scalar admits two different quantizations [21]. In 
this case, ip( x ' can either be a source or an expectation value according to the standard quantization or 
the alternative quantization, respectively. In our paper, we will only focus on the standard quantization. 
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The "shooting" method states that we can start with a initial value of ip at the tip 
r and then perform the numerical calculations of the EoMs of ip Eq.(19) provided that 
the infinite boundary conditions ip^ = are satisfied. At the critical point of the phase 
transition, the quantity of ip is very close to zero. Therefore, we have set the initial value 
of ip to be 0.01 in our numerical calculations. 
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Figure 2: The marginally stable curves of scalar fields corresponding to various critical 
chemical potentials in the cases of different squared mass. The critical chemical potentials 
for different curves are in the sequence /ij? luc < fi T c cA < fif een . 

Fig. (2) shows the multiple marginally stable curves of the scalar fields R(x) 3 in x = 1/r 
coordinates for various m 2 . Take the plot of m? = —15/4 for example, the first three 
lowest-lying chemical potentials are in the sequence /i)? lue < // c cd < /if cen . We find that 
the values of /i c 's obtained from the "shooting" method are perfectly consistent with the 
He's in Table (1) derived from the "marginally stable modes" method. The blue line is 
for the minimal value of \i c . It starts from a very small initial value at the tip x — 1 and 
then monotonically decrease to zero at the infinite boundary x — 0. There are no other 
intersecting points between the blue line and the R(x) = axis. This is the reason we call 
the mode with fi c rj 1.888 of node n = 0. For the red and green lines, they are related 
to fi c = 4.234 and \x c = 6.615 respectively. The red line has one intersecting point with 

3 Note that we have transformed the scalar field ip(r) to R(x) in x = 1/r coordinates. 
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R(x) = axis while the green line has two. Therefore, as we have mentioned above, we can 
call the two modes with nodes n = 1 and n = 2 respectively. However, the red and green 
lines are expected to be unstable, because radial oscillations in x-direction of R(x) will 
cost energy [18]. In addition, the arguments above are also appropriate for other diagrams 
in Fig. (2). 



3 P-wave perturbations 

In this section we extend our work to study an AdS soliton with p-wave vector fields. We 
consider a five- dimensional SU(2) Einstein- Yang-Mills theory with a negative cosmological 
constant following [22] . The action is 

S = J d*x^-g[\{R - A) - \F; u F a ""], (22) 

where F" is the field strength of the SU(2) gauge theory and F^ u = d^A^—d^A^+e^A^Al. 
a,b,c — (1, 2, 3) are the indices of the SU(2) Lie algebra generator. A ^ are the components 
of the mixed- valued gauge fields A = A^^dx 11 , where r a are the generators of the SU(2) 
Lie algebra with commutation relation [r a , r b ] = e abc r c . And e abc is a totally antisymmetric 
tensor with e 123 = +1. 

As a consistent solution of the system, in the probe limit, the background of the metric 
can also be an AdS soliton solution like (1) (We have scaled L = l,r = 1), 

dv 2 

ds 2 = ^r— - + r 2 (-dt 2 + dz 2 + dy 2 ) + r 2 g(r)dx 2 , (23) 
r 2 g(r) 

where we have set f(r) = r 2 g(r) = r 2 (l — 1/r 4 ). 
We adopt the ansatz for the gauge field as [18] 

A(t, r) = (f)(r)T 3 dt + i)(t,r)r 1 dz. (24) 

Note that in order to consider the QNMs of the vector field, we have assumed ip(t,r) 
depends on t and r. In this ansatz, the gauge boson with nonzero component ip(t, r) along 
z-direction is charged under A 3 = <f)(r). According to AdS/CFT dictionary, <p(r) is dual to 
the chemical potential in the boundary field theory while ip(t, r) is dual to the ^-component 
of some charged vector operator O. The condensation of ip(t,r) will spontaneously break 
the U(l) 3 gauge symmetry and induce the phenomena of superconducting on the boundary 
field theory. 

Define x = 1/r, the EoMs for <j)(x) and i/j(t,x) = e~ luJt R(x) in the x coordinate are 

0" + (- - ~W -—0 = 0, (25) 
g x g 

R" + ( 9 ---)R' + ^ ± ^R = 0. (26) 
9 x g 
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where a prime denotes the derivative with respect to x. A simple solution is that 4>{x) = 
Const. = fi and R(x) = which corresponds to an neutral AdS soliton with a constant 
electric potential. 

In order to find the marginally stable modes of R(x), we can follow the steps in the 
previous section. Eq.(26) can be transformed to 

S(x)d 2 x R(x) + I^-d x R{x) + JM-R( X ) = 0, (27) 
x — 1 (x — l) z 

where, the coefficient functions are 

S(x) = (28) 
x — 1 

T(x) = x 3 (-g + xd x g), (29) 

U(x) = x 4 (oo 2 + /2 2 )(x -1). (30) 

Polynomially expand these three coefficient functions to a finite order and R(x) to an 
infinite order, such as 

M 

S{x) = ^8 n (x-l) n , (31) 

n=0 

N 

R(x) = lim V a n (x- l) n . (32) 

N^>oo — ' 

71=0 

where M is a finite integer. Substitute (32) and (31) into Eq.(27) we reach a recursion 
relation as Eqs.(16) and (17) in the previous section. The solutions of u can be derived by 
requiring R(0) = 0. 

In the practical numerical calculations, we take the order of the expansion of R(x) to 
be N = 300. Notice again that in the potential U(x) (30) u and fi are symmetric. Using 
the trick we have adopted in the preceding section, we can easily find the marginally stable 
modes u = when /x = /i c . The first three lowest-lying critical chemical potentials are 

ti e » 2.265,4.742,7.156. (33) 

The minimal critical chemical potential /i c ~ 2.265 is in perfect agreement with the results 
in [10, 13]. 

In order to study the behavior of the field R(x), we should first know the boundary 
conditions for R(x) and then make use of the "shooting" method to numerically calculate 
it. Notice again that we still assume the electric field A t = Const. = /i here. At the tip, 
R(x) behaves as 

R( x ) = a + b{x - 1) + • • • , (34) 
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while 



R(x) = R {0) + R {1) x 2 + ■■■ . (35) 

at the infinite boundary. In the following calculation we will set = in order not to 
source the field theory on the boundary. 



0.010 

0.008 
0.006 
0.004 

R(x) 

0.002 
0.000 
-0.002 
-0.004 



Figure 3: The marginally stable curves of vector fields R(x) corresponding to various 
chemical potentials. The critical chemical potentials for different curves are in the sequence 

. blue ^ . . red ^ . green 
H'c ^ H'c ^ r'c 

We can see from Fig. (3) that there are multiple marginally stable curves corresponding 
to different critical chemical potentials. The blue line is related to the minimal value 
H c = 2.265 which is obtained from the "shooting" method. Once again, we find that this 
critical value of chemical potential is identical to the one derived from the "marginally 
stable modes" method. In addition, other values of fi c are also consistent. From what we 
have learned in the last section, the blue line corresponds to the modes of node n = 
while red and green lines are related to the modes of node n = 1 and n = 2 respectively. 
However, the last two modes, i.e., n = 1 and n = 2 are unstable due to the cost of energy. 

4 Discussions and Conclusions 

In this paper, we have studied the marginally stable modes for the s-wave and p-wave 
perturbations in AdS soliton background with a constant electric potential. At some critical 
chemical potentials, the marginally stable modes, i.e. u = will arise. The importance 
of the marginally stable modes is that they can reveal the instabilities of the background 
which in our paper means that the neutral AdS soliton will become unstable to develop 



P-wave 
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charged scalar (vector) "hairs" in this AdS soliton background. Although the detailed 
phase transitions cannot be seen in the study of QNMs, it actually has been announced 
in the previous works [6, 9, 10] by the study of thermodynamics such as the free energy. 
This phase transition in the gravity side will map to an insulator/superconductor phase 
transition on the boundary field theory. 

Despite that we do not exactly know the phase structures through the marginally stable 
modes, they can actually indicate the onset of the phase transition. This has been argued by 
Gubser in studying the holographic superconductor phase transitions [2, 18]. In particular, 
marginally stable modes can be obtained by studying the QNMs of the perturbations. The 
widely used method to study QNMs in asymptotically AdS spacetime was Horowitz and 
Hubeny's method [17]. In this paper, we have adopted this method to find that at some 
critical chemical potentials, there indeed appeared the marginally stable modes. These 
critical chemical potentials were in good agreement with those obtained from the previous 
studies [6, 10, 13]. We also took advantage of the "shooting" method to numerically 
plot the behaviors of the scalar (vector) fields in the radial direction. Therefore, one can 
intuitively see the "nodes" of the marginally stable modes. In addition to the modes of 
the minimal critical chemical potentials, we also studied other less lowest-lying marginally 
stable modes by both methods. We asserted that these less lowest-lying marginally stable 
modes are unstable because the oscillations of the scalar (vector) fields in radial direction 
will cost energy. 
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